Optically observable zero-age main-sequence 

O stars 

By NOLAN R. WALBORN 
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. A list of 50 optically observable O stars that are likely on or very near the ZAMS is presented. 

They have been selected on the basis of five distinct criteria, although some of them exhibit more 
than one. Three of the criteria are spectroscopic (He II A4686 absorption stronger than in normal 
luminosity class V spectra, abnormally broad or strong Balmer lines, weak UV wind profiles for 
their spectral types), one is environmental (association with dense, dusty nebular knots), and one 
is photometric (derived absolute magnitudes fainter than class V). Very few of these stars have 
■ been physically analyzed, and they have not been considered in the current framework of early 

massive stellar evolution. In particular, they may indicate that the earliest, embedded phases are 
not as large a fraction of the main-sequence lifetimes as is currently believed. Detailed analyses 
of these objects will likely prove essential to a complete understanding of the early evolution of 
massive stars. 
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1. Introduction 

It is often stated that zero-age main-sequence (ZAMS) O stars should not be and are 
not observed. This view arises from at least three sources: star-formation theory, which 
suggests that the embedded accretion (merger?) phases constitute a significant fraction 
' of the main-sequence lifetimes of massive stars (2.5 Myr for the most massive); statistical 

studies of UCHII and IR objects relative to optically observed ones; and detailed physical 



analyses of optical O-star samples that find very few on the ZAMS. For instance, Repo- 
lust, Puis, & Herrero (2004) analyzed 24 relatively bright O stars and found only one, 
HD 93128 in the Carina Nebula compact cluster Trumpler 14, on the ZAMS. However, 
selection effects may be contributing to this view. If the optically observable near-ZAMS 
phase of massive stars is relatively brief, it must be sought in very young regions, which 
k> , may be distant and/or extincted. Also, it is possible that some IR objects are no longer 

embedded, but rather viewed along unfavorable sightlines in galactic disks or through 
local, peripheral remanent dust clouds. For example, if we did not have such fortunate 
lines of sight toward NGC 3603 and 30 Doradus, we might be quite confused about their 
evolutionary status (Walborn 2002). Such must be the case for at least some objects. 

Over the past 35 years, this author and others have encountered numerous optical 
O stars that appear to be very young for various morphological reasons. These results 
are scattered throughout the literature and have not been generally recognized by star- 
formation and evolutionary specialists, or even by quantitative spectroscopists, whose 
analyses are essential for the former. It is hoped that this summary presentation will 
provide a useful stimulus toward rectifying the omission. It will be seen that many of the 
objects in question are in the Magellanic Clouds. But also, HD 93128 will reappear in 
the plot. 

2. Categories of candidate ZAMS O stars 

The current sample of 50 optically observable, likely ZAMS O stars is listed in Table 1, 
along with some normal standard stars for comparison. The ZAMS candidate list was 
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complete to the author's knowledge as of May 2006, although further candidates are being 
discovered as of this writing, e.g., in the SMC cluster NGC 346 by Evans et al. (2006), 
and in a new LMC fiber survey by I. Howarth (private communication). The sample is 
divided into five categories according to the distinct, principal discovery criteria, although 
many of them actually display more than one, as noted in the Comments. Three of these 
criteria are spectroscopic, one is environmental, and one is photometric. Discovery and/or 
data references are included in the Table. The five categories will now be discussed in 
turn. 

2.1. Vz spectra 

A luminosity classification for stars earlier than spectral type 09 was introduced by 
Walborn (1971, 1972, 1973a). It is based upon the selective emission effects (Walborn 
2001) in He n A4686 and the N in triplet AA4634-4640-4642, i.e., the Of phenomenon. 
These same lines display a negative effect in absorption with increasing luminosity in 
the MK 09-BO classification, which was hypothesized to be caused by filling in of the 
absorptions by the same emission effects producing the Of phenomenon, thus providing 
the basis for a luminosity classification at the earlier types. A luminosity sequence at 
spectral type 06.5 in modern digital data is shown in Figure 1. As can be seen there, 
this He II line is a strong absorption feature in class V, which then weakens, neutralizes, 
and finally comes into emission above the continuum in the la supergiant. (Correlatively, 
the N in, already weakly in emission at class V, increases in strength with increasing 
luminosity.) 

Walborn (1973) first noted that in the Trumplcr 14 O-dwarf spectra, the A4686 ab- 
sorption appeared stronger relative to the other He lines than in typical class V spectra. 
This very compact cluster in the Carina Nebula appears to be very young; it contains the 
02 If* prototype, likely pre-WN object HD 93129A (Walborn et al. 2002b). Penny et al. 
(1993) found that this cluster is approximately 550,000 years old, while Repolust et al. 
(2004) derived an age of 150,000 years for HD 93128. Subsequently, even more extreme 
examples were found in the LMC giant H n regions 30 Dor and Nil (Parker et al. 1992; 
Walborn & Parker 1992; Parker 1993; Walborn & Blades 1997). Some of the Nil spectra 
are reproduced in Figure 2, and some new Galactic examples kindly provided for this 
presentation by V. Niemela and N. Morrell are shown in Figure 3. 

A fairly obvious hypothesis is that the stronger A4686 absorption in these probably 
very young objects is an "inverse" Of effect, or more precisely, that typical class V spectra 
already have some emission filling in that line while these objects have less, and hence may 
be less luminous and less evolved, i.e., nearer to the ZAMS. To denote that hypothesis, 
the luminosity class notation Vz has been introduced. Some analytical support has been 
provided by the work of Venero, Cidale, & Ringuelet (2002), Martins et al. (2005), and 
Mokiem et al. (2006), but a homogeneous analysis of the full sample in high-resolution 
data with state-of-the-art photospheric/wind models is essential to investigate whether 
the Vz stars as a class have systematically higher gravities and lower luminosities than 
class V. In Table 1, the three blocks of Vz stars correspond to the Galaxy, LMC, and 
SMC. 

A misunderstanding of the Vz definition in late-0 spectra should be clarified here. 
At early O types, the He II A4686 absorption should be stronger than He II A4541. At 
type 07, He II A4541 is equal to He I A4471, so that A4686 is stronger than both in a 
Vz spectrum. At later O types, however, A4541 weakens more rapidly with advancing 
type than A4686 in normal class V spectra; thus, a late-0 Vz spectrum must have A4686 
stronger than He i XJ^ll. 
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He I I He II I He I I He II I He II 



4000 4200 4400 4600 

Wavelength (A) 

Figure 1. A luminosity sequence at spectral type 06.5. The rectified spectrograms are separated 
by 0.4 continuum units. The spectral lines identified below are He I AA4026, 4471 and He II 
AA4200, 4541, 4686. N III AA4634-40-42 emission is marked above HD 93146, likewise Si IV 
A4089 and S IV AA4486-4504 (Werner & Rauch 2001) in HD 163758. Note the comparable 
strengths of the AA4541, 4686 absorptions in the class V spectrum, and the weakening, then 
transition to emission of the latter with increasing luminosity, while the N ill emission increases 
smoothly. The Si IV absorption has a positive luminosity effect, which is more sensitive at later 
types. Courtesy of Ian Howarth. 

2.2. Vb spectra 

W. Morgan frequently remarked on a peculiarity of the Orion Nebula Cluster OB stars, 
namely broader Balmer lines than in normal class V spectra, with profiles that did not 
appear to be rotational; Morgan & Keenan (1973) reproduced photographic spectrograms 
of 9 1 Ori C and the 07 V standard 15 Mon to illustrate the effect. At some point, he 
introduced the notation Vb to denote such spectra, by analogy with lb for less luminous 
supergiants and Illb in Keenan's subdivision of late-type giants. (Abt 1979 cites a pho- 
tographic atlas prepared with Morgan that appeared in 1978 as the source of the new 
notation, but it does not appear there.) Abt (1979) and Levato & Malaroda (1981, 1982) 
presented further examples among B and A dwarf spectra in very young clusters. 

As for the Vz category, the question arises whether the Vb phenomenon might be 
caused by higher gravity and lower luminosity than in class V, since it is well known 
that the Balmer lines weaken with increasing luminosity in OB spectra due to the de- 
creasing Stark effect. A careful analysis of the five OB components of 9 1 and 9 2 Ori 
has been presented by Simon-Diaz et al. (2006), including a thorough investigation of 
line-broadening mechanisms and comparison with standard objects (including 15 Mon) 
analyzed with the same techniques. Interestingly, they find that H and He lines in the 



N. R. Walborn: Optically observable zero-age main- sequence stars 




Wavelength (A) 

Figure 2. Vz spectra in Lucke-Hodge 10/Henize Nil in the LMC (Walborn & Parker 1992). The 
ordinate ticks are separated by 0.25 continuum units. See the Fig. 1 caption for identifications of 
the He lines. Note the increased strength of the He II A4686 absorption relative to the other He 
lines in these spectra, and also the great strength of the Balmer lines in LH10— 3126. LH10— 3264 
is in the dense nebular knot N11A, which produces the very strong nebular emission lines. 
Courtesy of Joel Parker. 



Orion stars tend to be broader than in the best fitting models, most systematically in 
the BO. 5 V spectrum of 6 1 Ori D. However, they derive similar gravities to those of the 
comparison stars and find that the Orion stars are somewhat off the ZAMS, although 
uncertainties in the latter result remain because of the extinction law and, as they point 
out, the effect of initial rotational velocities on the location of the ZAMS. In Table 1, a 
value of R = 4.25 (average of 3.0 and 5.5) has been used to calculate the absolute mag- 
nitudes of the Trapezium stars, because the actual value in the Orion Nebula remains 
uncertain (Robberto et al. 2004); thus their subluminosity is uncertain as well. 

9 1 Ori C is now known as the first O-type magnetic oblique rotator (Donati et al. 
2002; Smith & Fullerton 2005; Gagne et al. 2005; Wade et al. 2006) and an extreme 
spectrum variable, including its abnormally weak UV wind profiles (Walborn & Nichols 
1994; Stahl et al. 1996). Simon-Diaz et al. (2006) provide a detailed discussion of the 
effects of the variability on the quantitative analysis. Moreover, there is now evidence 
that the spectrum of 15 Mon may be variable as well; Simon-Diaz et al. also discuss the 
effects of its spectroscopic companion in that connection. Further quantitative analyses 
incorporating these complications are indicated to ascertain the nature of the Vb Balmer 
profiles and the physical origin of the peculiarity. 
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Figure 3. New Galactic Vz spectra, rectified and separated by 0.2 continuum units. See the 
Fig. 1 caption for identifications of the stellar lines; in addition, N IV A4058 emission is marked 
in HD 64568 here. Courtesy of Virpi Niemela and Nidia Morrell. 



2.3. Weak winds 

It is now well known that the ultraviolet stellar-wind profiles in O-type spectra display 
strong correlations with the optical spectral types, including an increase in strength 
with increasing luminosity (Walborn, Nichols-Bohlin, & Panek 1985). That study also 
detected four main-sequence stars with abnormally weak wind profiles for their spectral 
types, including 9 1 Ori C already discussed in the previous section. HD 5005A in the 
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young cluster NGC 281 is also located in a trapezium system, while HD 42088 ionizes 
the small, dusty H n region NGC 2175, which might appear knot-like (cf. next section) 
if it were at the distance of the Magellanic Clouds. Again, the inverse behavior of these 
weak wind profiles with respect to normal spectra suggests that these stars may be less 
luminous and less evolved than typical class V objects, but quantitative confirmation 
of this hypothesis is so far lacking. Martins et al. (2005) found a relatively weak wind 
for HD 42088, but a normal HRD location; however, the distance of this star is highly 
uncertain. (It should be noted that most of the stars in Table 1 have not yet been observed 
in the far UV.) 

2.4. Nebular knots 

Walborn & Blades (1987) reported the detection of O-type spectra in two dense nebular 
knots within the 30 Doradus Nebula in the LMC, suggesting that they might correspond 
to very young objects just emerging from their natal cocoons. Subsequent work, both from 
the ground and with HST, has amply confirmed that suggestion and revealed additional 
examples in 30 Dor (Walborn & Blades 1997, Walborn et al. 1999, 2002a); several have 
been resolved into multiple (trapezium) stellar systems by HST. The very strong nebular 
emission lines in these objects obliterate the stellar He I absorption and thus preclude 
accurate classification, as only the He n can be seen, leading to the 03-6 V type (there 
is no uncertainty in the luminosity class, which depends on only He II A4686); in fact, 
several examples re-observed spectroscopically with HST, its very high spatial resolution 
suppressing the nebular background, turn out to have even later O types than that range 
previously assigned from the ground (Table 1). The reddening and absolute magnitudes 
of these stars have kindly been derived by L. Ubeda from the WFPC2 photometry using 
the CHORIZOS code (Mafz-Apellaniz 2004); a normal reddening law had to be adopted 
because of the limited wavelength coverage. 

LH10-3264 (Parker et al. 1992, Walborn & Parker 1992; Fig. 3 here), which is also 
Vz, is another interesting case in the LMC compact H n region Nil A; it has also been 
resolved into a multiple system with HST by Heydari-Malayeri et al. (2001), and the entry 
in Table 1 (N11A-7) corresponds to the brightest component. N81 is a similar object in 
the SMC, investigated with HST by Heydari-Malayeri et al. (2002); the weakness of the 
stellar- wind profiles in these stars may be caused by a combination of extreme youth and 
the SMC metal deficiency. 

2.5. Subluminosity 

In the spectroscopic study of the 30 Doradus stellar populations by Walborn & Blades 
(1997), the derived absolute magnitudes for a number of the O V stars fall below the 
calibration of that luminosity class (cf. their Fig. 3). Those that are not also classified 
as Vz are listed separately in Table 1. It is reasonable to hypothesize that they may 
be nearer to or on the ZAMS. However, a basic uncertainty, which also applies to some 
subluminous cases discussed in the previous sections, must be recognized: because of the 
lack of further information, the absolute magnitudes have to be derived with a normal 
value of R = 3. For a fixed distance modulus and variable extinction, a larger value of R 
would yield brighter absolute magnitudes. To resolve this issue, photometric observations 
covering a wider wavelength range (ideally, UV through IR) must be obtained to support 
the derivation of reddening laws toward each individual star. Indeed, there is evidence 
that the reddening law may vary among different lines of sight within a complex H n 
region environment, because of spatially diverse mixtures of dust grain properties. In any 
event, large values of R are associated with large dust grains and very young objects, so 
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Figure 4. Derived absolute visual magnitudes for candidate ZAMS and normal O stars plotted 
against their spectral types. The bold solid line is the ZAMS from Schaller et al. (1992), while 
the lighter solid lines are the corresponding isochrones for 1, 2, 3 Myr. The dashed line is the 
luminosity class V calibration from Walborn (1973a). The normal stars fall near the latter, while 
most of the ZAMS candidates lie below it. Courtesy of Leonardo Ubeda. 

these stars may be near the ZAMS in either case, but again, further observations and 
quantitative analyses are required to definitively establish their physical status. 

3. Discussion 

All of the absolute magnitudes from Table 1 are plotted against the spectral types 
in Figure 4, where they are also compared with the Schaller et al. (1992) ZAMS and 
isochrones, as well as the luminosity class V calibration (Walborn 1973). The different 
categories of ZAMS candidates are distinguished with different symbols. Subject to the 
uncertainties already discussed, the plot provides preliminary support for the present 
hypothesis: most of them fall below the class V calibration and near the ZAMS. Note that 
HD 93128 lies adjacent to the ZAMS, in agreement with the detailed analysis of Rcpolust 
et al. (2004). The normal class V stars cluster about the calibration, as well they should 
since they contributed to its derivation. Interestingly, the figure shows that a typical 
03 V star is 1 Myr old, 05 V 2 Myr, 06-7 V 3 Myr, and 08-9 V perhaps 4-5 Myr. Less 
luminous stars at a given type must be younger. It is also clear that unresolved multiple 
systems move points upward in the diagram (recall that O'.'l corresponds to 5000 AU at 
the LMC), while underestimated values of R move them downward. 

The present sample of candidate ZAMS stars is optimum for followup with high- 
resolution spectroscopy and state-of-the-art atmospheric/wind analysis. Beyond the ba- 
sic question of their ZAMS status or otherwise, it may be hoped that such study will 
elucidate any sequential relationships among the different categories, thereby advancing 
our detailed understanding of early massive stellar evolution. 
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